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We develop a theoretical framework to study slowly rotating compact stars in a rather general class 
of alternative theories of gravity, with the ultimate goal of investigating constraints on alternative 
theories from electromagnetic and gravitational- wave observations of compact stars. Our Lagrangian 
includes as special cases scalar-tensor theories (and indirectly f{R) theories) as well as models with 
a scalar field coupled to quadratic curvature invariants. As a first application of the formalism, we 
discuss (for the first time in the literature) compact stars in Einstein-Dilaton-Gauss-Bonnet gravity. 
We show that compact objects with central densities typical of neutron stars cannot exist for certain 
values of the coupling constants of the theory. In fact, the existence and stability of compact stars 
sets more stringent constraints on the theory than the existence of black hole solutions. This work 
is a first step in a program to systematically rule out (possibly using Bayesian model selection) 
theories that are incompatible with astrophysical observations of compact stars. 
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I. INTRODUCTION 

Compact stars as nuclear physics laboratories. 

Studies of compact stars in general relativity have been 
textbook material for decades P, Q ■ Neutron stars can 
be considered "cold" by nuclear physics standards, so 
their mass-radius relation M(R) is uniquely determined 
by the equation of state (EOS) of matter at high densi- 
ties, i.e. by the relation between pressure and energy den- 
sity P{p). From an observational point of view, one usu- 
ally assumes general relativity to be correct. Under this 
assumption (which of course is backed up by a wealth of 
observational evidence Q), the Holy Grail of astronomi- 
cal observations of neutron stars is the determination of 
the EOS from measurements of macroscopic properties 
such as masses, radii and moments of inertia. 

Better observational estimates of neutron star masses 
and radii are progressively improving our understanding 
of the EOS. Lindblom Q presented a concrete scheme for 
reconstructing P{p) from observations of M{R). More 
recently. Read et al. Q approximated the high-density 
EOS by piecewise polytropic models, showing that cur- 
rent astrophysical measurements yield stringent con- 
straints on the piecewise polytropic parameters. In the 
same spirit, Lindblom @ proposed to replace piecewise 
polytropes by spectral expansions, that should give a 
more faithful representation of the EOS. 

Our understanding of the functional form of the EOS 
from observed masses and radii has made imp ressive 
strides in the recent past_0-Q (see also [l^, lllj for re- 
views). Demorcst et al. |12l | recently determined a value 
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M = 1.97 ± 0.04 Mq for the mass of PSR J1614-2230, a 
pulsar in a white dwarf-neutron star binary system. This 
precisely measured mass is large enough to rule out many 
candidate EOSs [l^ [l^ . Vice versa, theoretical progress 
in microscopic calculations based on chiral effective field 
theory is leading to a better understanding of neutron- 
rich matter below nuclear densities, and hence to more 
stringent constraints on the mass-radius relationship [l^ . 

Compact stars as strong gravity laboratories. 

Most studies of the possibility of reconstructing the EOS 
from compact star observations assume that general rel- 
ativity is the correct theory of gravity. General relativ- 
ity passed all observation tests so far but the "real" 
theory of gravity may well differ significantly from it in 
strong field regions. In fact, cosmological observations 
and conceptual difficulties in quantizing Einstein's theory 
suggest that general relativity may require modifications. 

Compact stars are an ideal natural laboratory to look 
for possible modifications of Einstein's theory and their 
observational signatures [l6| . Besides ruling out specific 
models for the EOS, experiments may (and should) try 
to rule out also alternative theories of gravity that are 
unable to explain observations. A comprehensive study 
of how EOS models and alternative theories affect macro- 
scopic observable quantities of compact stars requires a 
Bayesian model selection framework, where one compares 
the predictions of any specific theory of gravity (and of 
different EOS models) against the growing body of ob- 
servational data. Of course, an important prerequisite 
of any such analysis is the construction of stellar models 
in the largest possible family of alternative theories of 
gravity that are not ruled out by weak-field experiments, 
cosmological constraints or observations of compact bi- 
nary systems. The present work is a first step in this 
direction. 
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The plan of the paper is as follows. In Sec. |TT1 to 
put our work in context, we briefly review some stud- 
ies of compact stars in alternative theories of gravity. 
In Sec. mil we present the Lagrangian for what we call 
"extended scalar-tensor theories" . In Sec. IIVI we write 
down the equations describing the structure of static and 
slowly-rotating stars in this class of theories, and we dis- 
cuss our chosen models for the EOS. In Scc.|V]we present 
numerical results and discuss their implications. We con- 
clude by discussing extensions of the present work to 
other theories and comparisons with observations. Unless 
stated otherwise, we use geometrical units {G = c = 1). 



II. STARS IN ALTERNATIVE THEORIES: A 
BRIEF REVIEW 



The study of compact stars in alternative theories of 
gravity has a long history. In this paper we begin a sys- 
tematic exploration of stellar structure in a large class of 
modified gravity theories. Far from providing a compre- 
hensive review, here we point to some relevant literature, 
mainly to put our work in context. 

i) Scalar-tensor theories. Scalar-tensor gravity is one 
of the simplest and best-motivated modifications of gen- 
eral relativity, because scalar fields are predicted by al- 
most all attempts to incorporate gravity into the stan- 
dard model [17[ • Therefore it should come as no surprise 
that most work on stellar structure concerns variants of 
scalar-tensor theory. The equations of hydrostatic equi- 
librium in the best-known variant of scalar-tensor theo- 
ries (Brans-Dicke theory) were first studied by Salmona 
(isj . Soon after, Nutku [l^ explored the radial stabil- 
ity of stellar models using a post-Newtonian treatment. 
Hillebrandt and Heintzmann [20j analyzed incompress- 
ible (constant density) configurations. Zaglauer [2l| car- 
ried out a detailed calculation of the so-called "sensitiv- 
ities" of neutron stars, which determine the amount of 
dipolar gravitational radiation emitted by compact bina- 
ries in scalar-tensor theories [22j. Most of these studies 
found that corrections to neutron star structure are sup- 
pressed by a factor 1/wbd, where wbd is the Brans-Dicke 
coupling constant. At present, the most stringent bound 
on this parameter (wbd > 40, 000) comes from Cassini 
measurements of the Shapiro time delay Q. 

As pointed out by Damour and Esposito-Farese j23j . 
the coupling of the scalar with matter can produce a 
"spontaneous scalarization" phenomenon by which cer- 
tain "generalized" scalar-tensor theories may pass all 
weak-field tests, and at the same time introduce macro- 
scopically (and obscrvationally) significant modifications 
to the structure of compact stars. More detailed stud- 
ies of stellar structure [U, [25| , numerical simulations of 
collapse j26l - l28j and a stability analysis [1^ confirmed 
that "spontaneously scalarized" configurations would in- 
deed be the end-state of stellar collapse in these the- 
ories. In fact, spontaneously scalarized configurations 
may arise as a result of semiclassical vacuum instabilities 



[3 Of . Tsuchida et al. [3l| extended the Buchdahl inequal- 
ity {Al / R < 4/9 for incompressible stars) to generalized 
scalar-tensor theories. For a comprehensive study of an- 
alytic solutions and an extensive bibliography, see [s^ . 
ii) f(R) theories. Theories that replace the Ricci scalar 
i? by a generic function f{R) in the Einstcin-Hilbert ac- 
tiorl^ can always, at least in principle, be mapped into 
scalar-tensor theories [35|. The existence of compact 
stars in metric f{R) models that have been proposed to 
explain cosmological observations, such as the Starobin- 
sky model [s^ , was studied by many authors [37l - l4l] | with 
controversial results. One possible explanation of the 
partial disagreement between different authors is that the 
mapping between f{R) theories and scalar-tensor theo- 
ries is in general multivalued, and therefore one should be 
careful when considering the scalar-tensor "equivalent" of 
an f{R) theory [l^. A perturbative approach to stellar 
structure in f {R) gravity is also possible [43| . 
Hi) Higher- curvature gravity. Besides theories where 
the Lagrangian is a generic function of R, it is of interest 
to consider theories where the Lagrangian is built out of 
quadratic [i^l (or even higher-order) contractions of the 
Riemann and Ricci tensors. As we explain below, the 
requirement that the field equations should be second- 
order means that quadratic corrections must appear in 
the Gauss-Bonnet (GB) combination 
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where Rabcd is the Riemann tensor and Rab is the Ricci 
tensor. Since the GB term in four dimensions is a topo- 
logical invariant, the GB combination introduces mod- 
ifications to general relativity only when coupled to a 
nonzero scalar field or other forms of matter. The sim- 
plest and better motivated cas§j is Einstein-Dilaton- 
Gauss-Bonnet (EDGB) gravity [i^l, where the GB term 
is coupled to a dynamical scalar field, the dilaton. The 
EDGB correction to the Einstein-Hilbert action appears 
in low-energy, tree-level effective string theory |46| . 

The study of EDGB gravity in relativistic astrophysics 
has been limited to a mathematical analysis of black hole 
solutions fisl. lirl - ligj and, more recently, to their possible 
observational signatures [H, [s^, HH . To our knowledge, 
the present study is the first investigation of compact 
stars in the theory. Static black holes in EDGB gravity 
only exist when a > H^. Hence, we shall restrict our 
study to the case of positive a. 

iv) Parity-violating theories. Chcrn-Simons gravity 
is the simplest theory that allows for parity- violating cor- 
rections to general relativity [Hlj. Due to the nature of 



Here and below we refer to f{R) theories in the metric formal- 
ism. Theories of the Palatini type have conceptual problems: for 
example, spherically symmetric polytropic "stars" present cur- 
vature singularities [331. l34l| . 

In analogy with f{R) models, /(T^-qb) models have been studied 
in a cosmological context. Observational constraints on 
models are quite tight (see e.g. Sections 12.3 and 12.4 of 1351 ) 
and we will not include them in our analysis. 
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the Chern-Simons corrections, all spherically symmetric 
solutions of Einstein's theory are also solutions of Chern- 
Simons gr avity . However, spinning objects in the non- 
dynamical |53| and dynamical [s^l versions of the theory 
are affected by the Chern-Simons coupling. Future ob- 
servations of the moment of inertia of compact stars may 
strongly constrain the parameters of the theory [U • 

v) Lorentz-violating theories. Einstein-aether theory 
introduces a dynamical unit timelike vector coupled to 
gravity as a natural way to implement Lorentz violation 
in Einstein's theory. In the parameter space compati- 
ble with Solar System constraints, spherically symmetric 
neutron stars in Einstein-aether theory have a lower max- 
imum mass than in general relativity [ssl [56| . Another 
popular Lorentz-violating theory is Hofava gravity. The 
matching conditions necessary to obtain stellar solutions 
in this theory were considered in (57| , but (to our knowl- 
edge) there are no phenomenological studies of compact 
stars using realistic EOS models. 

vi) Massive gravity. Recently, Damour et al. |58| re- 
considered the discontinuity problem of massive gravity 
and its possible resolution through Vainshtein's nonlinear 
resummation of nonlinear effects. As part of this study, 
the authors investigated the viability of spherically sym- 
metric stars in the theory. They showed that some solu- 
tions show physical singularities, but also that there exist 
regular solutions interpolating between a modified gen- 
eral relativistic interior and a de Sitter exterior, with cur- 
vature proportional to the square of the putative graviton 
mass. A more phenomenological study of observational 
constraints (including stellar rotation) is still lacking. 

vii) Eddington inspired gravity. Banados and Fer- 
reira [59j recently proposed a theory that is equivalent 
to general relativity in vacuum, but differs from it in the 
coupling with matter. An interesting aspect of this the- 
ory is that singularities cannot form in early cosmology 
and during gravitational collapse jsol [60j . The maximum 
mass of compact stars in the observationally viable sec- 
tor of Eddington-inspired gravity may be larger than in 
general relativity, even for "ordinary" EOS models [60t . 

viii) Gravitational aether, f(T), TeVeS and other 
theories. Some alternatives to general relativity that 
were proposed to explain cosmological observations have 
also been analyzed, at least to some extent, in the 
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context of compact stars. Among these theories we 
can list "gravitational- aether" theory [6lj, /(T) gravity 
[6^ and Bekenstein's TeVeS [g^l- In higher-dimensional 
braneworld models, the embedding of four-dimensional 
stellar solutions "on the brane" within acceptable higher- 
dimensional solution is a nontrivial problem [6^ (but see 
[gsI for related work in a slightly different context). 



III. EXTENDED SCALAR-TENSOR THEORIES 

From the previous summary it should be clear that it is 
nearly impossible to discuss all strong-field modifications 
of general relativity in a unified framework. However, in 
this section we show that on the basis of some rather gen- 
eral arguments we can easily write down a Lagrangian 
encompassing the first four classes of theories reviewed 
above (namely general scalar-tensor theories, f{R) theo- 
ries, EDGE gravity and Chern-Simons gravity). 

Our starting point is a Lagrangian in which gravity 
is coupled to a single (generically charged) scalar field 
(j) in all possible ways, including all linearly independent 
quadratic curvature corrections to general relativity. We 
call these models "extended scalar-tensor theories" . The 
most general Lagrangian of such a theory contains several 
functions of the scalar field in the combination 

c= hM)R--imda(i^*d''cj,-vM) 

+ /l(|0|)i?2 + hm)RabR'''' + hmRabcdR'^'"'''' 

+ /4(|0|)i?a6c/i?'''^'' + /:mat [*,A2(|0|)g,,] , (2) 

where is the dual of the Riemann tensor, which in- 

troduces possible parity- violating corrections j52| . From 
the Lagrangian above, the equations of motion read: 



Gat+T" \Hab 
JO 



- Jab + K^ab] 



1 

Wo 



(0) 
b 

(3) 



where T^*"* = 2(-.g) ^/'^6Sm/Sgab is the matter stress- 
energy tensor in the Jordan frame, 

T^if = 7 [29(a0*96)0 - gabdccl>*d'cf\ - gabV 

+2VaV6/o - 2.g,,V2/o , (4) 
and. following the notation of Ref. [3], we have defined 



-Hab = -^v[l^Vb)R - 2RV^avi]^ + gab \2RV'vi'^ + 4t-(i) V^i?j + h 



2RabR ' 2VabR - ^gab [R^ - ^OR] 



(2) 



lab = -t'gVb)i? - 2t.(2) [V^aRbf - V^i?a6] + VS'^)^^^ _ 2i?,(,VS, 
^ ./2 [2i?^'*i?acfc!i — Vafci? + ORab + ^.9ah [OR — RcdR'^''']] , 

'.ab] + mUV'vf - h [2 [RabR - 4.R'^Racbd 



9ab 



Jab = -8u^^' [V(ai?b)'' - VRabj i -^-ii-acb ^ 

JCab = H'^e-deia'^^Rb)" + 4Vrfi;W*i?(/b)' , 



VabR - 2nRab] - \gab [R^ - ^RcdR''']] 



with v^a* = Va/i(|0|), Vah = V^Vb, □ = V^V" and e'''"^'' the Levi-Civita tensor. The modified Klein-Gordon 
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TABLE I. Specific models obtained from tlie Lagrangian ([2}. Here k = {16nG) 



fo fi f2 fs fi 



V 



1 



A 



General relativity k 
Scalar-tensor (Jordan frame) [24] F((\)) 
Scalar-tensor (Einstein frame) [2J| k 
/(7?) [35] 

Quadratic gravity [44] 
EDGE [45| 
Dynamical Chern-Simons [54] k 
Boson stars [66] k 
























K 



as,<l> a4(t> 

fi 

I3(j> 

ij 











1 

2k 
2k 

1 
1 
1 
1 



1 
1 

A{4>) 



fo 



■1/2 



1/2 



perfect fluid 
perfect fluid 
perfect fluid 
perfect fluid 

perfect fluid 
perfect fluid 
perfect fluid 




equation reads 

— Js^abcd^ — J4^abcd n — A A 1 J , (t)J 

together with its complex conjugate. In the equations 
above, a prime denotes a derivative with respect to 

As shown in Table HI this theory is sufficient to discuss 
stellar structure in many of the alternative theories listed 
in Sec. HI] (and it can also describe boson stars in general 
relativity, if we work in vacuum). As a matter of fact, 
some terms in the Lagrangian ^ are redundant. For 
example, in ordinary scalar-tensor theories the functions 
/o and 7 can be removed via a conformal transformations 
of the metric and a redefinition of the scalar field; i.e., 
by reformulating the theory in the Einstein frame [23 ]. 
However, depending on the explicit form of /o and 7, 
these transformations can be hard (if not impossible) to 
write in a closed analytic form. For this reason we find 
it convenient to start from the general Lagrangian ([2]), 
which reduces to standard scalar-tensor theories in the 
Jordan frame if A( ](/)]) = 1 and in the Einstein frame if 
/o(l</)l) = l/(167r) and 7(10]) = 1. Another advantage of 
this approach is that, at least in principle, it should also 
encompass generic f{R) theories, which are equivalent to 
particular scalar-tensor theories (but see [42j for possible 
issues with this point of view). 



A. Simplifying the model 

For generic coupling functions, the terms V-abi lab, Jab 
appearing on the left-hand side of the equations of mo- 
tion introduce higher-order derivatives of the metric func- 
tions, unless quadratic terms in the curvature enter the 
action in the GB combination (HI). The GB combination 
corresponds to setting /2 = — 4/i and = fi in our 
model. Thus, if we only want second-order equations of 



motion the Lagrangian ^ must reduce to 

£= /o(101)i? + /l(|</>|)7^^B + hm)Rabcd*R'''''' 

-7(l0l)5a0*9"<A - Vm + £„,at A\\c^\)gab\ .(7) 
In order to avoid the complications related to higher- 
order derivatives, from now on we will specialize to this 
Lagrangian. 

IV. PERFECT FLUID COMPACT STARS IN 
EXTENDED SCALAR-TENSOR THEORIES 

A. Static solutions 

We begin by looking for static, spherically symmetric 
equilibrium solutions of the field equations with metric 

dsl = -B(r)dt^ + — — + r'^dO^ + sin"^ 9 dip^ 

1 — 2m(r)/r 

and a charged, spherically symmetric scalar field 

^{t,r)=<f{r)e-^-K (8) 

Because of the assumed spherical symmetry, the Pontrya- 
gin density vanishes {Rabcd* R"'^'^'^ = 0) and the equations 
of motion do not depend on f^. Our ansatz for the scalar 
field also implies that the Klein-Gordon equation ([5]) and 
its conjugate coincide. 

We consider perfect-fluid stars with energy density p{r) 
and pressure P{r) such that 

rZt ^ r;,-:,,,,,,,, = (p + P) + g^^'P , (9) 

where the fluid four-velocity = (1/a/B, 0, 0, 0). Note 
that the matter fields are defined in the Jordan frame. 
The stress-energy tensor in the Einstein frame reads 

T^iv' = J^^<l>\)Tiiv To close the system, as usual, we 
must also specify an EOS P = P(p). 

The field equations for a static, spherically symmetric 
perfect-fluid star in extended scalar-tensor theories read 
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£^00 = TT^ r'V($) -r^^($)2p + 4rVo(«')m' -4r''/o($)$' + 6r^7i/o($)$' + 16rm/{($)$' 

B 

-48mV{(*)*' + 2r'*/^($)m'$' - 16r^f[{<i>)m'^' + 48rm/{ ($)m'$' - r^7($)$'^ + 2r*mj{^)^'^ - 2r^$'Vo (^) 
+4r^m$'Vo (*) - 16r2m$'V{'(*) + 32rm2$'V"($) - 2r(r - 2m) [rVo(*) + 8m/{($)] $" = , (10) 
£;ii = r*cj27($)$2 _ (r - 2m)S' {2r2/o($) + [?'^/o(«') - 8(r - 3m)/{($)] $'} + rB {4/o($)m + r [r'^A{^)'^P 

-r2y($) + (r - 2m)$' (-4/^($) + r7($)$')] } = , (11) 
Scons = 4r^A{^)B^fo{'P)PA'{^)^' + r''A{^fB [/o($) ((P + p)B' + 2SP') - 2BP/^($)$'] + $' x 

{r/o(<i>)(r - 2m)B'^ P/^($) + 8m/{(<&)] + 2B {r/^($) [-r4c^2^(<i>)<i>2 + (r - 2m)P' (rVoCf) 
-8(r - 3m)/{($)) $'] + /o($) [-r3mP7^($) + 8rmP7{($) - 2Am^B'f[{^) + r"^ B' f[^{<^>)m' 
-8r^B'f[{^)m' + 2ArmB' J[{^)m' + r^uj^^'^H^) + r^-f{^) (2ra;2$ + (r - 2m)P'$') 
-r(r - 2m) (rVo(*) + 8m/{($)) B"] } + 2r2p2 {r/^($) [r'^V{<^) + (r - 2m)$' (4/^($) - r7($)$')] 
+/o($) [r (4/^($)m' + r {-rV'{<^) + r7'($)$'2 + 27($) (- (-2 + m') + r$"))] - 2m (2/|^($) 
+r (r7'($)$'2 + 7($) (3$' + 2r$"))))] } = , (12) 
-Escai = r-(r - 2m)S'2 [rVo($) + 8m/{($)] + 25 {B' [r'^/^($) (3m + r (-2 + ?n')) + 8(r - 3m)/{($) (m - rm')] 
-r^w2$2^'($) + r'*7($) [2ra;2^ + (r - 2m,)P'$'] - r(r - 2?7i) P/o(*) + 8m/{($)] B") 

+2r^B^ {r2A($)3(3P - p)A'($) + 4/^($)m' - rV($) + [-27($) (3m + r (-2 + ?7i')) + r(r - 2m)7'($)$'] 
+2r(r-2m)7($)$"} = 0, (13) 



where Pqo and En arc the {0, 0} and {1,1} components 
of the modified Einstein equations, i?cond = ^ aT""^ = 
and i?scai denotes the field equation for the scalar field. 

To construct spherically symmetric and static stellar 
configurations, we must solve the system above imposing 
regularity conditions at the center of the star, i.e. 

m(0) = 0, p(0) = pc, $(0) = $c, $'(0) = 0. 

(14) 

More in general, any field can be expanded close to the 
center as 

X{r) ^ + X(iV + X(2)r2 + ©(r^) , (15) 

where X schematically denotes any of the variables p, 
P, ^, B and m. By using the field equations, all coef- 
ficients X*-*' ultimately depend on two parameters only, 
say = pc and ~ (f>c- Finally, the value of $c is 
fixed through a shooting method in order to obtain an 
asymptotically flat solutioi:0: $ — ?> as r — )• oo. 

The outcome of this shooting method is a one- 
parameter family of solutions characterized only by the 
central density pc- For any value of pc, we can compute 
the mass M and the radius Rs of the star. As usual, the 
mass is obtained from the asymptotic behavior at infinity 

2M 

P(r) -> 1 , (16) 



This condition can be easily generalized to allow for nonvanishing 
values of the scalar field ($ — > $oo 7^ as r — > 00), like those 
arising in f{R) theories and in other modified gravity theories 
that want to reproduce cosmological dynamics. 



whereas the radius is computed by imposing the usual 
matching condition at the stellar surface, P{Rs) = 0. In 
the exterior of the star, P = p = 0. Finally, the baryonic 
mass, corresponding to the energy that the system would 
have if all baryons were dispersed to infinity, is defined 
as 

m = m,jd^.^gu^nir), (17) 

where n denotes the baryonic number. The normalized 
binding energy fh/M — 1 is positive for bound (but not 
necessarily stable) configurations. 

B. Slowly rotating models 

Once a static stellar model is known, it is easy to con- 
struct the corresponding slowly rotating model by gener- 
alizing the classic work by Hartle [g^I . For this purpose, 
we consider the metric ansatz 

ds^ = -B{r)dt^ + [1 - 2m(r)/r]"^ dr^ + r^dO^ 
+ sin^ B {dLp - ^dtj 

^ dsl - C{ry sm"^ Odtdifi + OiC'^) . (18) 

The stress-energy tensor for a rotating fluid can be easily 
constructed from Eq. © and the four-velocity 

= {u*,0,0,l]u*} , = ^-{gtt + 2nt^ + n^g^^), 

(19) 

where 51 is the angular velocity of the fluid. 
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Note that, for any nonconstant /4(|(^|), the gravito- 
magnetic part gttp of the metric would source scalar per- 
turbations through first order parity-violating terms in 
Eq. dni). If /o,/i,A and 7 are constant and F = 0, 
the field equations admit the same spherically symmet- 
ric solutions as in general relativity, and the background 
scalar field vanishes. In this case, the only first-order 
corrections arise from the {t, Lp\ component of the Ein- 
stein equations and from the (perturbed) scalar equation, 
since the stress-energy tensor is quadratic in the scalar 
field. These two equations can be solved for C,{r) and 
for the scalar field perturbation for slow rotation and 
small scalar fields. This was done in Ref. [HJl, where the 
authors studied slowly rotating neutron stars in dynam- 
ical Chern-Simons gravity (/4(|(^|) = al^D- When (as 
in our case) the background scalar field is nonvanishing, 
all nontrivial Einstein equations would acquire first-order 
corrections. Hence, if Ad^l) is not constant, Hartle's 
formalism cannot be applied. In particular, other first- 
order metric corrections (in addition to the gravitomag- 
netic part gtip) must be included in Eq. (|18p . This would 
result in a system of equations which is very difficult to 
decouple. To avoid these issues, from now on we will 
simply assume that /idt/f)!) =const. 

Using the definition of the four-velocity ([T^ and lin- 
earizing in the angular velocity $7, we find that the solu- 
tion to the field equations corresponds to the background 
(nonrotating) solution plus the solution C(r) of an ordi- 
nary differential equation coming from the {t, Lp\ compo- 
nent of the Einstein equations, namely: 

C"(r)+Ci(r)C'(r) = [C(r)-r!] X (20) 
(r - 2A/) [r2/o($) - 4(r - 2A/)/( ($)$'] ' 

where 



2rB{r - 2M) (r2/o($) - 4(r - 2A/)/{($)$') 
[r(r - 2M)B' (4(r - 2M)/{($)$' - r^/ol*)) 
+2B (-r2/o($) {7M + r {M' - 4)) + (r - 2Af)x 
+ X ($' (rVo(*) + 12/{(*) [M + r (A/' - 1)) 
-4r(r - 2Af )$'/('($)) - 4r(r - 2Af )/{($)$"))] . 

Equation (j20p must be solved by imposing regularity con- 
ditions at the center of the star: C(0) = Cc, C'(0) — 0. We 
must also require continuity of C(r) at the stellar radius. 
The asymptotic behavior at infinity reads 

2 / 

C ^ Coo + ^ , (21) 

where J denotes the angular momentum. For the solu- 
tion to be asymptotically flat, we must impose Coo = 0. 
This can be easily achieved by noting that Eq. (PH)) is 
invariant under the transformation 

C^C-^, n^n~T^, (22) 



where rj is some constant. Therefore we can proceed as 
follows [g^]: (1) integrate Eq. (^0)) imposing regularity at 
the center and extract Coo and C,'^^ at some large (but fi- 
nite) radius Too] (2) find the physical value of the angular 
velocity, i.e., $7 — Coo- After this translation, C, -> 2J/r^ 
at infinity; (3) compute 

J = _ lim . (23) 

r— >oo 6 

As we vary 17 — Cc we obtain models with different 
specific angular momentum J/Af^. As long as 17 <C 
\J AI/R^, the slow-rotation approximation is consistent. 
Ignoring terms of 0(17^), the moment of inertia is given 
by / = J/fl and it does not depend on Cc, but only on 
the stellar mass. Therefore we need to integrate Eq. ([20)1 
only once in order to obtain / for a given mass. 

With the slowly rotating solution at hand, we can also 
study the possibility of ergoregion formation [g^ ■ The er- 
goregion can be found by computing the surface at which 
gtt vanishes, i.e., from Eq. (fT8|): 

- B{r) + C^r^ sin^ e = . (24) 

On the equatorial plane we simply have 

rC(r) = /B(0 , (25) 

and, due to the linearity of the field equations, C will 
scale linearly with 17. Thus, one needs only a single in- 
tegration in order to compute the zeros of Eq. ([?5|) as 
functions of 17. For a given value of 17, there can be no 
zeros (i.e. no ergoregion), two distinct zeros (with the er- 
goregion located between them) or two coincident zeros. 
The "critical frequency" at which we have two coinci- 
dent zeros, say 17c, is the minimum rotation frequency 
for which an ergoregion exists. The slow-rotation ap- 
proximation imposes 17 < 17,ns, where the mass shedding 
frequency is defined as 17„is = y^M/R^, following Har- 
tle's conventions. The existence of an ergoregion requires 
17 > 17c, so an ergoregion can exist (in the slow-rotation 
approximation) only if 

17c < ^ms ■ (26) 

C. Equation of state 

We wish to establish some limits on the parameters 
of allowed theories of gravity, given a set of EOS mod- 
els compatible with our present knowledge of nuclear 
physics. Instead of general relativity being assumed, the 
parameters characterizing any specific theory will be con- 
strained based on astrophysical observations. These con- 
straints will be sensitive to our assumptions on the EOS, 
that we describe in this section. 

A list of the EOS models used in this work is given 
in Table [Hi For code testing purposes, we have con- 
sidered the same polytropic model used by Damour and 
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TABLE II. EOS models used in this work. 



EOS 


Reference 


Polytropic 


[23] 


FPS 


[69] 


APR 


[70] 


Causal limit 


[15, 7^ 



Esposito-Farese [23 



p = nmii + K 



nomb I n 



r - 1 \no 



P — KnQjni) 



' g, no = 0.1 fm 



no, 

(27) 
2.34 and 



with nil, = 1.66 x 10 '^'^ ■ 
K = 0.0195. 

We also considered two nuclear-physics motivated 
models (FPS jl^ and APR [t^], in the standard nomen- 
clature), which are respectively a soft EOS and a more 
standard realistic EOS, as well as for the stiffest possi- 
ble EOS constructed by combining the upper limit in the 
crust-core transition re gion of Hebelcr et al. 



causal limit EOS as in 71 



151 with a 



The polytropic model (P7)) 
gives results which are quantitatively very similar to 
those for FPS EOS. 

We must remark that the FPS EOS seems to be ruled 
out by the recent observation of a neutron star with 
M = (1.97 ± 0.04) Mq [ill, at least if we limit consid- 
eration to nonrotating model^ within general relativ- 
ity. However, these observations could be explained in 
terms of modified gravity at high density, rather than 
by invoking a different EOS. In fact, in some alterna- 
tive theories the maximum mass of a neutron star can 
be sensibly larger than in general relativity [13, [sO] ■ An- 
other important motivation to use the FPS EOS is to 
make direct comparison with previous work. We explic- 
itly checked that our two independent codes (written in 
Mathematica and C-| — h) are in excellent agreement 
with Refs. [^, [t^ [tJ] in the general relativistic limit[3. 



V. COMPACT STARS IN GAUSS-BONNET 
GRAVITY 

As a first application of the formalism discussed above, 
in the remainder of this paper we study neutron stars in 
EDGB gravity. We defer a more general study of the full 
theory derived from the Lagrangian ([7]) to future work. 



* Rapidly rotating neutron stars have a larger maximum mass 
which, using the FPS EOS, is still marginally compatible with 
the observational errors in a small region of the parameter space: 
some rapidly rotating models in Table 3 of 1721 have a gravita- 
tional mass compatible with the value measured in [T3 |. 

^ Incidentally, the moment of inertia shown in Fig. 3 of l54i is 
not computed using the FPS EOS, as erroneously written in the 
caption of that figure. This explains why our results for the FPS 
EOS do not agree with those in lEi . 



EDGB gravity is obtained from the Lagrangian ([7]) by 
considering a real scalar field = $ (or uj = 0), fo = k = 
(167r)-\ I/ = and 



/i 



167r 



(28) 



where a and /3 are coupling constants. Static black holes 
in EDGB gravity only exist when a is positive j45j and 
for this reason we will only consider a > 0. When /3 = 
•^2, this theory arises as a low-energy correction to the 
tree-level action in heterotic string theory [i^. Here we 
adopt a phenomenological point of view and consider a 
and (3 as free (real) parameters. We will show by an 
explicit calculation that, under reasonable assumptions 
for the nuclear EOS, the observation of compact stars 
with certain observed properties (such as mass, radius 
or moment of inertia) leads to the existence of rather 
stringent exclusion regions in the two-dimensional (a , /3) 
parameter space. 

Some results are shown in Figs. [T][3] for different EOS 
models and different values of a and /3. In each figure, 
the top two panels show the mass-density relation and the 
mass-radius relation for static (nonrotating) stars. The 
bottom-left panel shows the binding energy as a function 
of the central density. In the bottom-right panel we dis- 
play the moment of inertia as a function of (gravitational) 
mass. 

In our numerical calculations, we observed that the 
scalar field in the interior of the star is always small: in 
special cases it can be as large as $ '--^ 10~^, but more 
typically $ < 10~^ in most of the parameter space. In 
the small-field limit, $ ^ 1, the coupling /i in Eq. 
can be Taylor-expanded: 



167r/i($) 



(29) 



responding to a = 20M|,, jS'^ = 1 and a = 10M|, = 4 

20M|, and indeed 



Since the first term is constant and the GB term is a topo- 
logical invariant, the first nonvanishing corrections arise 
from the second term. Therefore, in the small-field limit 
the equilibrium structure depends only on the product 
af5 of the coupling constants. This is confirmed by our 
numerical results in Figs. [T][3l for instance the lines cor- 
20M2 

both correspond to the same a/3 
they lie almost exactly on top of each other. A similar 
degeneracy will occur for any other functional form of 
fi{^), provided that the scalar field remains small ev- 
erywhere, so that a Taylor expansion similar to Eq. ([25)1 
holds. In this sense, most of our results remain valid for 
a generic function /i($); and not only for EDGB gravity. 

It is clear from Figs. [l]|3]that, regardless of the EOS, 
for a > the coupling to the dilaton tends to reduce 
the importance of relativistic effects. Indeed, as shown 
in Fig. m the maximum gravitational mass Mmax mono- 
tonically decreases as a function of the product a/3 of the 
EDGB coupling parameters. Thus in EDGB gravity (as 
well as in general relativity) soft EOS models, like FPS, 
should be ruled out by observations of high-mass neutron 
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FIG. 1. Compact star models in EDGB gravity for different values of the parameters a and /3, using the APR EOS. In the 
bottom-right panel we show the recent observation of a neutron star with M « 2Mq and a possible future observation of the 
moment of inertia confirming general relativity within 10% [t^. Curves terminate when the condition (|30p is not fulfilled (cf. 
also the exclusion plot in Fig. [Sj . 



stars. This is similar to what happens in gravitational- 
aether theory [6l| and in Einstein-aether theory [56| . 

For small values of the product a/3, the maximum mass 
in Fig. 21 corresponds to a local maximum in the mass- 
density relation (cf. the upper left panels of Figs. [T]|31). 
In general relativity these local maxima (or, equivalently, 
inversion points in the mass-radius diagram) correspond 
to marginally stable equilibrium configurations, and so- 
lutions to the right of the first maximum are unstable to 
radial perturbations (see e.g. (H). We conjecture that the 
same property should hold also for extended scalar-tensor 
theories. This was proved for particular self- grav itating 
configurations involving scalar fields [2^, [tI, l77j , but a 
more detailed stability analysis would be desirable (see 
also the discussion in [HI). In EDGB theory, spheri- 
cally symmetric solutions can be constructed only up to 
a maximum central density pJf'^(a,/3), for reasons ex- 
plained below: see in particular the discussion around 
Eq. (|30p . For large a/3 this maximum central density is 
such that the first local maximum in the mass-density 
curve is never reached. In Fig. HI all values to the left 
of the solid circles correspond to a maximum mass ob- 



tained from the radial stability criterion. Values to the 
right of the solid circles correspond instead to the mass 
Mniax obtained at the critical value of pc beyond which 
we cannot find spherically symmetric perfect fluid solu- 
tions anymore. 

A. Constraints on the EDGB couplings 

In the near future, observations of the double pulsar 
may provide measurements of the moment of inertia to an 
accuracy of ~ 10% [zl] (but see Ref. [ill for some crit- 
icism). Furthermore, precise observations of the mass- 
radius relation could be obtained from thermonuclear X- 
ray burst [t^, [sO] . These observations could be used in 
the context of a Bayesian model-selection framework to 
place strong constraints on EDGB gravity and, more gen- 
erally, to remove the degeneracy between different EOS 
models and different proposed modifications of general 
relativity. 

Nevertheless, even without assuming any particular 
EOS, we can set rather stringent theoretical constraints 
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FIG. 2. Compact star models in EDGB gravity for different values of the parameters a and /?, using the FPS EOS. Curves 
terminate when the condition (|30fl is not fulfilled (cf. also the exclusion plot in Fig. [5]). 



on the EDGB parameters. Indeed, as shown in Figs.[T][31 
depending on a and /3, there is a maximum value of the 
central density pc^ above which no compact star models 
can be constructed. For a given central density, the criti- 
cal value of a/3 can be computed analytically in the small 
$ limit, as follows. We first compute the series expan- 
sion ([15]) up to O(r^). The resulting expressions are not 
very illuminating, but in general the series coefficients 
contain square roots, whose argument must be positive 
to ensure the existence of physical (real- valued) solutions. 
When $c ^ 1, by imposing this "reality condition" wc 
find 



1 



77767rP4pc 



(30) 



For a given value of a/?, the condition above implies that 
a maximum central density, p^^^ ^ exists. 

Equation ([50)1 is in good agreement with numerical re- 
sults, as shown in Fig.O This figure is basically an exclu- 
sion plot: it shows the maximum allowed values of a/3 as 
a function of the maximum central density p™^" for differ- 
ent EOS models and nonrotating stars. For small values 



of a/3 (i.e., on the right of the figure), a local maximum 
in the mass-density relation is reached and the maximum 
central density p^^^ corresponds to the local maximum of 
M(pc); i-e. to the first inversion point in the mass-radius 
relation. If our stability conjecture is correct, no sta- 
ble static configurations can be constructed in the region 
above these lines. The local maximum is never reached 
to the left of the points marked by filled circles. In this 
case, the M{pc) curves terminate before reaching a local 
maximum, and the maximum central density p™'*'^ is sim- 
ply the point where the cquilbirum sequence terminates. 
Dotted lines correspond to the analytical prediction (|30p , 
which agrees very well with the numerical value at which 
we cannot compute equilibrium models anymore. The 
bottom line is that no static models (either stable or un- 
stable) can be constructed in the shadowed regions above 
these exclusion lines. 

The quadratic EDGB corrections are expected to be 
stronger in high-density (high-curvature) regions, so the 
most stringent bounds should come from the stiffest EOS 
models. Indeed, among the models we consider, the 
strongest and weakest constraints come from the Causal 
EOS and from the FPS EOS, respectively. 

The bounds on the central density can be translated 
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Pc (lO^Vcm^) M/Mq 

FIG. 3. Compact star models in EDGB gravity for different values of the parameters a and /3, using a causal EOS. Curves 
terminate when the condition (|30[l is not fulfilled (cf. also the exclusion plot in Fig. [5]). 



into constraints on the maximum mass, which is an ob- 
servable quantity. As shown in Fig. |4l for a given EOS 
the maximum mass is a monotonically decreasing func- 
tion of a/3. 

The requirement that the maximum mass A/max sup- 
ported by the theory should be larger than some trusted 
observed value, can place a direct upper bound on a/3. In 
Table Uni we consider /W^ax > IAMq, A/^ax > 1-7Mq 
and Afn^ax ^ 1.93Mq (which is the lower bound of the 
recent observation in [l2|) and we translate them into 
upper bounds on a/3 using the data plotted in Fig. |4l 
The Causal EOS models predict an unrealistically large 
maximum mass, so all neutron stars with Mmax ^ 2.8Mq 
would place very mild constraints on alternative theories, 
and we omit this EOS from Table Hill 



TABLE III. Constraints on the EDGB parameters from an 
observation of a nonrotating neutron star with mass M. For 
the values of M we consider, constraints using the Causal 
EOS would allow values of a/3 larger than WOMq. 

EOS Afn,ax > 1.4Af0 A/^ax > 1.7A/q Afn,ax > l.QSM© 



FPS < 30.1Af| 
APR Q/3 < 50.3M| 



a/3 < 13.9M|, 
a/3 < 41.9Mq 



no models 
a/3 < 33.6M| 



For the "standard" value of the EDGB coupling (/3 = 
V2), if we consider the APR EOS as our "best candidate" 
EOS for a (nonexotic) neutron star interior within gen- 
eral relativity, the results in Table Hill imply a < 23.8/Vfg. 

This bound should be compared to the bound on a that 
comes from requiring the existence of black hole solutions 
in the theory [i^ [Hoj- For /3 = \/2, this requirement 
implies 



Ml 



< 70 



AI 



BH 



(31) 



where Mbh is the black hole mass. The observation of 
black holes with Mbh ~ 8Mq (such as Cyg XI) con- 
strains a < AAMq. The constraints on a coming from 
observations of compact stars (cf. Tabic HlI)) are already 
smaller by a factor ^ 2 than those coming from the exis- 
tence of stellar black holes, and they could become even 
more stringent in the near future. 

We may hope that future observations of the moment 
of inertia could place even tighter bounds on the theory. 
Unfortunately this seems unlikely. To understand why, 
we can either look at the bottom right panel of Fig. [1] 
or at Fig. [51 where we show the moment of inertia for 
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aPIMl 

FIG. 4. Maximum mass as a function of the product a/3 of the 
EDGB couphng parameters, for different EOS models and in 
the nonrotating case (cf. the main text for details). To the 
left of the filled circle, this maximum mass corresponds to the 
radial stability criterion; to the right, it corresponds to the 
maximum central density for which we can construct static 
equilibrium models. The recent measurement of a neutron 
star with M « 2Mq is marked by a horizontal line. Only 
the combination a/3 is bounded, due to the approximation 
$ < 1: cf. Eq. (Uni). 



for EDGB gravity, the most stringent constraints on a/? 
should come from mass measurements, rather than from 
measurements of the moment of inertia. 




FIG. 6. Moment of inertia normalized by its value in general 
relativity for the APR EOS at fixed values of the gravitational 
mass. Curves terminate at the bounds listed in Table[TTI] (cor- 
responding to the filled circles). The deviations from general 
relativity are always smaller than 5%. 




FIG. 5. Exclusion plot in the (p™^", a/3) plane for nonrotat- 
ing models. No compact star solutions can be constructed in 
the region above the dotted lines (cf. Eq. pop ). In the region 
above the thick lines (labeled by "RI", for "radial instabil- 
ity"), static configurations should be unstable against radial 
perturbations: see the main text for details. Filled circles in- 
dicate the maximum central density of radially stable stars in 
general relativity. 



the APR EOS (normalized by its value in general rela- 
tivity) as a function of a/3 for fixed values of the stellar 
mass. As it turns out, for values of a/3 smaller than 
those listed in Table IIIIl the moment of inertia can de- 
viate from the general relativistic value by 5% at most. 
The precision of future observations is expected to be 
~ 10% in optimistic scenarios [tHI- Therefore, at least 



Let us mention, for completeness, that we also studied 
the possibility of the formation of an ergoregion for slowly 
rotating stars in EDGB gravity. For our 'realistic" EOS 
models the condition ((26|) is never met, and therefore no 
ergoregion can form outside the star. This is because, for 
phenomenologically viable parameters, the relativistic ef- 
fects in this particular theory are actually smaller than in 
general relativity. We can anticipate that other sectors of 
the general theory described by the Lagrangian ([7]) could 
enhance relativistic effects and favor the existence of the 
ergoregion, with important implication for the stability 
of these solutions [6g. A more detailed analysis will be 
presented elsewhere. 



VI. CONCLUSIONS AND OUTLOOK 

Neutron stars are very promising laboratories to con- 
strain strong-curvature corrections to general relativ- 
ity. New proposed theories of gravity are usually tested 
against weak field observations and cosmological data, or 
by studying the existence and nature of black hole solu- 
tions. Our main goal in this paper was to develop a for- 
malism for a comprehensive study of stellar structure in a 
broad class of alternatives to Einstein's general relativity. 
We focused on a class of theories ( "extended scalar-tensor 
theories") where quadratic curvature corrections, non- 
minimal couplings and parity- violating terms are coupled 
to standard gravity through a single scalar field. Partic- 
ular cases of this model include, but are not limited to, 
quadratic gravity, EDGB gravity, generic scalar-tensor 
theories and f{R) theories (via their correspondence with 
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scalar-tensor theories). We wrote down the field equa- 
tions for static and spherically symmetric perfect-fluid 
stars in the general well as the leading-order cor- 

rections in a slow-rotation expansion. For a given model 
and a given central density, the formalism allows us to 
obtain the mass, radius, binding energy and moment of 
inertia of compact stars. In future work we will show 
how these theoretical predictions can be compared to ob- 
servations in order to constrain the parameter space of 
"extended scalar-tensor" and other alternative theories. 

As a first application of the formalism we studied stel- 
lar structure in EDGB gravity. We found that, in gen- 
eral, the GB coupling tends to reduce relativistic effects in 
compact stars. We also showed that there is an exclusion 
region in the two-dimensional plane of the GB coupling 
parameters beyond which no compact star solutions can 
be constructed: cf. Eq. ((SO]), Fig. |4] and Fig. [H 

Stability requirements for static models and future ob- 
servational data could constrain the theory even further. 
The existence of high-mass neutron starss put the most 
stringent constraints on EDGB gravity (cf. Table IIIip . 
As it turns out, these bound are tighter (by a factor of a 



few) than the bound coming from the existence of black 
hole solutions in EDGB theory, given in Eq. (jSlJ. They 
are also tighter than the bounds that could come from 
future precision measurements of the moment of inertia. 

In this sense, to our knowledge, the existence of large- 
mass neutron stars provides the best constraint on the 
EDGB coupling parameters obtained so far. Further ex- 
plorations of stellar structure and better observational 
data on the mass-radius relation (see e.g. [zl, HO]) have 
the potential to exclude a larger region of the parame- 
ter space of alternative theories. There is of course the 
possibility that theoretical and observational work may 
give us hints on how to modify general relativity to make 
it compatible with the standard model, which would be 
even more exciting. 
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